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ABSTRACT 

We investigate the radial colour gradient and two-dimensional (2D) colour properties of 
22 E+A galaxies with 5.5A < H(5 EW < 8.5A and 49 normal early-type galaxies as a control 
sample at a redshift of < 0.2 in the Second Data Release of the Sloan Digital Sky Survey. We 
found that a substantial number of E+A galaxies exhibit positive slopes of radial colour gradi- 
ent (bluer gradients toward the centre) which are seldom seen in normal early-type galaxies. 
We found irregular 'Colour Morphologies' - asymmetrical and clumpy patterns - at the cen- 
tre of g~r and r—i 2D colourmaps of E+A galaxies with positive slopes of colour gradient. 
Kolomogorov-Smirnov two-sample tests show that g— r and r—i colour gradient distributions 
of E+A galaxies differ from those of early-type galaxies with a more than 99.99% signif- 
icance level. We also found a tight correlation between radial colour gradients and colours, 
and between radial colour gradients and 4000A break in the E+A sample; E+A galaxies which 
exhibit bluer colour or weaker I34ooo tend to have positive slopes of radial colour gradient. 
We compared the GISSEL model and E+A's observational quantities, US EW, I?4ooo and 
u—g colour, and found that almost all our E+A galaxies are situated along a single evolution 
track. Therefore, these results are interpreted as E+A galaxies evolving from H(5 EW ^ SA to 
H(5 EW ^ 5 A, with colour gradients changing from positive to negative, and with the irregu- 
lar 2D colourmap becoming smoother, during a time scale of ~300 Myr. Our results favor the 
hypothesis that E+A galaxies are post-starburst galaxies caused by merger/interaction, having 
undergone a centralized violent starburst. 

Key words: galaxies: general 



X 



1 INTRODUCTION 

IPressler & GunnI il983l Il992h discovered galaxies with myste- 
rious spectra while investigating high redshift cluster galaxies. 
These galaxies had strong Balmer absorption lines with no emis- 
sion in [Oil]. They were named "E+A" galaxies, since their 
spectra resembled a superposition of those of elliptical galax- 
ies (Mg5i75, Fe527o and Ca3934,3468 absorption lines) and those 
of A-type stars (Strong Balmer absorption)^. Since the life- 
time of an A-type star is about 1 Gyr, the existence of strong 
Balmer absorption lines shows that these galaxies have experi- 
enced starburst within the last Gyr. However, they show no sign 
of on-going star formation as non-detection in the [Oil] emis- 
sion line indicates. Therefore E+A galaxies are interpreted as post- 
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^ Because the spectra of elliptical galaxies are characterize d by K stars, 
these galaxies are sometimes called "K+A" galaxies (e.g., lFranxlll993t 
iDressler et alll999l:lBartholom ew. Rose & Gaba 200 1|. Following the first 
discovery, we refer to them as "E+A" throughout this work. 



starburst galaxies, that is, galaxies which 
starburst I Dressier & Gunn 1983, 1992; 


lave undergone truncated 
Couch & Sharpies 19871 


MacLaren. Ellis & Couch 1988.: 


Newberrv. Boroson & Kirshnen 


1990l:lFabricant. McClintock & BautJ 1 99 lU Abraham et alJl996l). 



At first, "E+A" galaxie s were found in cluster regions 



both in low redshift clusters jFranxlll993l ICaldwell et al 



Caldwell & Ros^ Il997t ICastander et alj l200lL iRose et al. 

and high redshift clusters ( Sharoles et al. 1985; Lavery <S 



199; 



2001) 



& Henr^ 

1986; Couch & Sharpies 1987; Broadhurst, Ellis & Shanfa 
1988; Fabricant. McClintock & Bautz 1991; Belloni et al 
1995: .Bargeret aL.1996: .Fisher et al...l99&: .Morris et ah. 199a 
ICouch et alJ Il998t TpTessler et alJ Il999h . Therefore, a cluster 
specific phenomenon (e.g.. Goto et al. 2003b) was thought to 
be responsible for the violent star formation history of E+A 
galax ie s. A ram-pre s sure st ripping model jS mtzsr&Baadj 
195ll : 'Gunn & Got? '197^; iFarouki & Shacird Il98d; lKent| 
19811: lA badi. Moore & Bower 1999; Fuiita & Nagashima 199^ 
Ouilis. Moore & Bowe r 2000; Fuiita 2004; Fuiita & Goto 2004) as 
well as tides from the cluster potential (e.g.. Fuiita 2004) may first 
accelerate star formation of cluster galaxies and later turn it off. 
However, recent large surveys of the nearby universe found many 
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E+A galax ies in the field regions jGotdl2003i Eooi iGoto et all 
l2003cl: IOuintero 20o4- At the very least, it is clear that these E+A 
galaxies in the field regions cannot be explained by a physical 
mechanism that works in the cluster region. E+A galaxies have 
often been thought to be transitional objects during cluster galaxy 
evolu tion, evolving phen omena such as the Butcher- Oemler effect 
fe.g.. iGoto et alJ 12003 j) . the morphology-density relation (e.g., 
iGoto et alJl2003el ). and the correlation between various properties 
of the galaxies with the environment (e.g., Tanaka et al. 20041). 
However, explaining cluster galaxy evolution using E+A galaxies 
may not be realistic anymore. 

One possible explanation for E+A phenomena is dust- 
enshrouded star formation, where E+A galaxies are actually star- 
forming, but emission lines are invisible in opti cal waveleng t hs due 
to heavy obscuration by dust. As a variant, Pog gianti & WiJ J2000I) 
presented the selective dust extinction hypothesis, where dust ex- 
tinction is dependent on stellar age, since the youngest stars in- 
habit very dusty star-forming HIT regions while older stars have 
had time to migrate out of such dusty regions. If O, B-type stars 
in E+A galaxies are embedded in dusty regions and only A-type 
stars have long enough lifetimes ('^l Gyr) to move out from such 
regions, this scenario can naturally explain the E+A phenomena. 
A straightforward test for these scenarios is to observe in radio 
wavelengths where the dust obscuration is negligible. At 20cm ra- 
dio wavelengths, the synchrotron radiation from electrons acceler- 
ated by supemovae can be observed. Therefore, in the absence of 
a radio-loud active nucleus, the radio flux of a star-forming galaxy 
can be used to estimate its curren t massive s tar formatio n rate (SFR) 
|Condoa.l992. : .Kennicutt.l99& iHopkins et al. 2005) . ISmail et alj 

1 19991) performed such a radio observation and found that among 8 
galaxies detected in radio, 5 have strong B aimer absorption with no 
detection in [Oil]. They co ncluded that massive stars are currently 
forming in these 5 galaxies. |Owen et ai J IT999I) investigated the ra- 
dio properties of galaxies in a rich cluster at 2; ~ 0.25 (A2125) 
and found that optical line luminosities (e.g., Ha+[NII]) were of- 
ten weaker than one would expect for the SFRs implied by the ra- 
dio emission. Miller & Owen ( 2001) observed radio continua of 15 
E+A gala xies and dete cted moderate levels of star formation in only 

2 of them. iGotoliioM^ performed 20cm radio continuum observa- 
tion of 36 E+A galaxies, and none of them are detected in 20cm, 
suggesting that E+A galaxies are not dusty-starburst galaxies. 

Alternatively, galaxy-galaxy interaction has b een known 
t o trigger star formation in the pair of galaxies jSchweizen 

!98^;'Laverv & Henrv* 1988^, 'Liu & Kennicutt^995a''b; SchweizeJ 
996; Nikolic, Cullen & Alexander 2004). Oegerle. Hill & Hoessel 
1991) found a nearby E+A galaxy with a tidal tail feature. High 
resolution Hubble Space Telescope imaging supported the galaxy- 
galaxy interaction scenario by showing some post-starburst (E+A) 
galaxies in high redshift clusters as having disturbed or inter- 
ac^ng signatures (Couch etal. 1994, 1998; Dressier et al. 1994; 
lOemler. Dressier & Butcheill997l) . lLiu & Kennicutll (1995.3 b) ob- 
served 40 merging/interacting syst ems and found that some of their 
spectra resemble E+A galaxies. iBekki. Shiova & CouchI i200lh 
modeled galaxy-galaxy mergers with dust extinction, confirming 
that such syst ems can produce spe ctra which evolve into E+A spec- 
tra. Recently, looto et alj i2003d) ; lOoto (2005) found that young 
E+A galaxies have more companion galaxies within 100 kpc, pro- 
viding strong support for the merger/interaction origin of E+A 
galaxies. 

As is mentioned above, most of the previous work is focused 
on the global properties of E+A galaxies. Needless to say, investi- 
gation of the internal properties of E+A galaxies such as age and 




Figure 1. SDSS r-band image of one of the nearest E+A galaxies 
(left;2=0.034) captured from the SDSS DR2 sample. The image size is 86.3 
arcsec X 66.9 arcsec. Note the dramatic tidal tails. 



metallicity gradients is also indispensable in solving the E+A mys- 
tery. Do E+A galaxies leave traces of merger/interaction? Where 
do the traces remain, in a centralized or decentralized starburst? 
The elucidation of these questions is essential in testing E+A evo- 
lution scenarios. Norton et al. (2001) performed long-slit spectro- 
scopic observation of 21 E+A galaxies, and found young stellar 
populations of E+A galaxies are more centrally concentrated than 
older populations, and old co mponents of E+A galaxies conform 
to the Faber-Jackson relation. 'Bartholomew . Rose & Gab^ i200lh 
also reported that E+A galaxies on average tend to have slightly 
bluer ra dial gradients towa rd the centre than do normal early-type 
galaxies. lYang et all J2004h presented HST observations of the five 
bluest E+A galaxies with « ~ 0.1 and reported details of disturbed 
morphologies and detected compact sources associated with E+A 
which are consistent with the brightest clusters in nearby starburst 
galaxies. Although these are important results, previous work often 
lacked statistical significance, since E+A galaxies are extremely 
rare. In addition, it is important to understand the spatial proper- 
ties of E+A galaxies along an evolutionary sequence. Furthermore, 
there is a cont amination problem in the E+A sample selected with- 
out using Ha. iKennicuttI il992ald) shows, Hq is the best star for- 
mation indicator in the optical wavelength since it is a strong line 
and it has fewer uncertainties (e.g., dust extinction, self-absorption, 
metallicity dependence) than the other lines. Indeed, iGotol i2003) 
reported that a n E+A sample wit hout Ha cut-off contains ~50% 
contamination. Isiake et alJi2004h also found that the criterion us- 
ing H(5 and [Oil] leads to a significant sub-population of disk sys- 
tems with detectable Ha emissions. We stress that a clean sample 
selection is essential in order to study a rare population of galaxies 
such as E+As. 

In this paper, we use publicly available true E+A galaxies 
(without Ha nor [Oil] emission) selected from the Sloan Digital 
Sky Survey (SDSS, York et al. 2000; Early Data Release, Stoughton 
et al. 2002, hereafter EDR; First Data Release, Abazajian et al. 
2003, hereafter DRI; Second Data Release, Abazajian et al. 2004, 
hereafter DR2) by Goto ( 2005). Both broad-band imaging and the 
spectroscopic survey of 10,000 deg^ of SDSS provides us with the 
first opportunity to study E+A galaxies in a very large number. Gotol 
( 2005) analyzed ~250,000 galaxy spectra in the DR2, and the num- 
ber of E+A homogeneous galaxies reached 266. The SDSS imag- 
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ing is somewhat poor (typical seeing size is ^--^1.5") compared with 
8m-class telescopes or HST, but this sample includes very nearby or 
large E+A galaxies like FigureQ which exhibits dramatic tidal tails 
or conspicuous disturbed morphologies. We then select 22 E+A 
galaxies with z < 0.2, large apparent size, and more strict criteria 
5.5A < HS EW, -l.OA < Ha EW and -2.5A < Oil EW (a 
positive sign is absorption), and investigate spatial properties; 2D 
colourmaps and g—r and r—i radial colour gradients. The radial 
colour gradients are compared with spectroscopic properties, and 
we compare these properties with evolution scenarios using SED 
models in order to understand the evolutionary sequence of E+A 
galaxies. 

This paper is organized as follows: In section |2| the defini- 
tions of two galaxy samples are summarized, as well as the spec- 
troscopic and imaging data. In section |3| we explain the reduction 
of the imaging data, e.g. an image convolution, position matching 
between passbands, a review of the concentration index and a def- 
inition of radial colour gradients. In section|4| we briefly describe 
our E+A morphologies, and study 2D g—r and r—i colour prop- 
erties, 2D colourmaps and radial colour gradients, with relations to 
other photometric/spectroscopic properties. We then investigate the 
evolution scenario for E+A radial colour gradients using a compar- 
ison between SED models and our observational data. In section 
|5| we discuss our results with other E+A-related studies. Our con- 
clusions are presented in section |6| The cosmological parameters 
adopted throughout this paper are Ho — 75kms~^Mpc~'^, and 
(n^, Ha, fife) = (0.3,0.7,0.0). 



2 DATA SAMPLES 

Our E+A galaxi es are selecte d from the publicly available cata- 
log described in iGotd i2005h . The mother sample of this cata- 
log contains ~250,000 objects classified as galaxie s with spec- 
trosco pic information in SDSS DR2 ( Abazaiian et al fcoO^ . TGotd 
llOOS) selected 266 E+A galaxies as those with 5.0A < H5 EW, 
-3.0A < Ha EW and -2.5A < [Oil] EW (absorption lines 
have a positive sign) after redshift (> 0.01) and S/N (>10 per 
pixel) cut-off. 

We have restricted our targets to those galaxies which have 
a redshift z of less than 0.2 and an apparent size of aeo (Pet- 
rosian 60% flux semi-major axis; described in section|3} more than 
2.82" to maintain enough sampling points for radial colour gradi- 
ents. This selection leaves 50 E+A galaxies. To minimize contam- 
ination of our E+A galaxies with remaining star formation activ- 
ity, we set more strict spectroscopic criteria 5. 5 A < H(S EW and 
-l.OA < Ha EW, and finally, we are left with 22 E+A galaxies 
of greater purity. 

We also randomly selected 50 normal early-type galaxies 
within z < 0.2 as a control sample from objects classified as 
galaxies in the SDSS DR2 catalog. The normal early-type galax- 
ies are selected using the elliptical-based concentration index, Ce 
(Yamauchi et al. 2005), a ratio of Petrosian 50% flux semi-major 
axis to 90% flux semi-major axis (See section|3); Ce < 0.33 for 
early-type galaxies ( Shimasaku et al. 2001). Although the standard 
(inverse) concentration index C = rso /rgo defined with the Pet- 
rosian flux in the circular apertures is significantly affected by in- 
clination, the e ffect of inclination is r emoved with the use of ellip- 
tical apertures ( Yamauchi e t al. 2005). The cut-off size of aeo is the 
same as that of E+A galaxies. But this subsample contains 1 galaxy 
with a deblending problem, which we eliminated. Finally, we use 
49 normal early-type galaxies as a control sample. 



The spectroscopic data, H(5, Hq and [Oil] equivalent width 
(EWs) and their errors a re measured by the flux-summing method 
described in lGotoH2005h (See also Goto et al. 2003c). The values of 
4000Abreak(D40oo) in our paper are the reciproca l numbers of val- 
ues in the SDSS catalog. We use DR2 atlas images I Stoughto n et alJ 
l2002h as imaging data to calculate 2D colourmaps and radial 
colour gradients of g—r and r—i, and concentration index Ce - The 
photometric system, imaging hardware and astrometric calibration 
of SDSS are des c ribed i n detail elsewhere ( Fukugita et al. 199^ 
Gunn etaljll993 ; iHogg et alJl200lL ISmith et alJl2002i; IPieret^ 
■2003,V 



3 DATA REDUCTIONS 

From the SDSS database, we can easily download the deblended 
galaxy images, called 'atlas images', by fpAtlas file in which im- 
ages are archived using Rice compression. In the SDSS images, 
however, the seeing sizes differ from passband to passband, and 
the relation between physical and imaging positions of a passband 
on an Atlas image is not equal to others, but slightly differs. There- 
fore we cannot directly compose g- and r-band images, or r- and 
i-band images using the coordinates given by the fpAtlas file. We 
adjust the positions of the multi-band images in the following way: 
First, we smear all images except those with the worst seeing 
(Jmax among the passbands to equalize the seeing sizes of all pass- 
bands. Assuming that the point-spread function is Gaussian, images 
are convolved with a Gaussian function of ctc, 

CTc = \/ Crgiax - (1) 

where a is the actual seeing of each image derived from the 
psf .width parameter in tsField files: 

a = (psf_width/1.06)/2.35, (2) 

whe re psf _width is equivalent to 1.06 FWHM for a Gaussian pro- 
file <Fan et all2003l) . 

Next, each passband image is resampled using positional in- 
formation, "cole" and "rowc"(real) in tsObj file, and "offsets 
wrt^ reference colour"(dcol and drow; integer) and "bounding 
box"(cmin and rmin; integer) in "master mask" in fpAtlas file (see 
source code of readAtlasImages). Where "cole" and "rowc" are 
at the exact centre of the object in fpC coordinates, and (cmin + 
dcol, rmin + drow) is the fiducial point of atlas images in fpC co- 
ordinates. The typical errors of "cole" and "rowc" are 0.01 ~ 0.1 
pixel(0.004 ~ 0.04 arcsec). 

In the case of making g—r 2D colourmaps without resampling 
r-band images, the preliminary position offset (xOpg,yOpg) of g- 
band pixels needed for Ix g-band resampling in an atlas image 
coordinate is 

^Op9 ~ (colcg — eolOg) — (cole,. — eolOr) 
yOpg — (rowcg — rowOg) — (rowe,. — rowOr) , 

where colO = cmin + deol and rowO — rmin + drow. And we 
consider the values {dOx, dOy): 

dOx = XOpg - [xOpg] 

dOy = yOpg - [yOpg] 

where [x] is the floor function and {dOx,dQy) indicates the distance 
from pixel grid and the resampled position of (;-band pixels (That 
is, resampling is not required when {dOx, dOy) = (0, 0)). 

^ with reference to 
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Resampling only g-band images is not reasonable. Both r- and 
p-bands should be resampled and optimized. Then we determine 
the position offsets by dQx, 



xQog = xQpg + (-d0:,/2) 

xQor = {-dOx/2) 

xOog = xOpg + (1 - d0^)/2 

xOor = (l-d0,)/2 



or 



and by dOy, 

yOog = 
yOor = 



yOo 



yOpg + {-dOy/2) 

{-dOy/2) 



yOpg + (1 - d0y)/2 

(l-d0,)/2 



for dOx 0.5 (5) 



for 0.5 < dOx , (6) 



for dOy <= 0.5 (7) 



for 0.5 < dOy . (8) 



These position offsets minimize declines in the resolution of r- and 
g-band images. The position offsets of u-, i- and z-bands, xQou, 
xQoi and xQoz, are also derived by equation jsj or j6j, and 
or (SJ, replacing g with u, i or z. The xOpu, xOpi and xOpz are 
calculated by equation ^3), replacing <; with u, i or z (Note that 
equation J4} is applied for the g-band only). In the case of r—i 
2D colourmaps, xloi, yloi, xlor and ylor are calculated using the 
same convention. When performing 1 x resampling, bilinear filter- 
ing is adopted. 

Then we compute the surface brightness SB of each pixel on 
a 1 X resampled image a nd correct the reddening due to dust ex tinc- 
tion in our Galaxy, using Schlesel, Finkbeiner & Davis ( 1998). The 
SDSS magnitudes m are derived by the following equations used 
in the SDSS photometric pipeline (PHOTO; Lupton et al. 2001), 



0.4(Uaa+I'kk-"air) 



'^asinh 



'^Pogso 



2.5 
In 10 



10' 



ln& 



-2.5 ■ log(///0) 



(9) 



(10) 



(11) 



where A'^cnt is counts within an aperture on the atlas image, Tcxp 
is exposure time in the fpC file, Uaa , Vkk and Uair are zeropoint, 
extinction coefficient and airmass in the tsField file, and b is the 
asinh softening parameters. See the DR2 Photometric Flux Cali- 
bration page jhttp://www.s dss.org/dr2/algorithms/fluxcal.html I for 
details of computing magnitude from SDSS CCD images. Finally, 
we calculate 'pixel by pixel' _ftr-corrected surface brightness us- 
ing the kcorrect.v3_2 library (Blanton et al. 2003). The same K- 
correction procedure is applied to each pixel in our galaxy based 
on the u,g,r,i and z colours of the pixel. The if-correction code 
uses several templates obtained as eigen vectors from the real SDSS 
spectra. During the fitting procedure, the code finds the bestfit out 
of the liner combinations of these template spectra. Therefore, the 
code does not try to break the age-metallicity degeneracy, but uses 
typical spectra with typical ages and metallicities to find the rest- 
frame colour of galaxies. Details of the A'-correction code itself are 
written in Blanton et al. ( 2003). 

The concentration index Ce is computed using an r-band im- 
age using elliptical apertures which eliminates the effect of the 
apparent axis ratio of the galaxy. The details of the effect and 
m ethod of deriving the a xis ratio and position angle are described 
in I Yamauchi et alji2005h . To calculate the concentration index Ce 
for the ellipse, we consider area Ae (a) of the ellipse of the semi- 
major axis a and axis ratio a, and the integrated flux Fe (a) within 



Ae{a). The Petrosian semi-major axis ap for a given 77 is defined 
by 

_ {fe(1.25ap)-Fe(0.8ap)}/{Ae(1.25ap)-Ae(0.8ap)} 

F4ap)/Ae{ap) ' 

where we take rj = 0.2, and the elliptical Petrosian flux Fp as 

Fp = F4kap) (13) 

with k set equal to 2, following the SDSS defiirition iStrauss et alJ 
l2002h . The Petrosian half-, 60%- and 90%-light semi-major axes 
150, 160 and ago are defined in such a way that the flux in the ellip- 
tical apertures of these semi-major axes are 50%, 60% and 90% of 
the elliptical Petrosian flux: 

F,(a5o) = 0.5Fp, Fe{am,) = 0.6Fp, F^ago) = 0.9Fp . (14) 



We define our concentration index Ce by 



(15) 



The axis ratio and the position angle used above are also adopted 
for the calculation of the radial colour gradient. 

The radial colour gradients are based on our pixel-to-pixel K- 
corrected g—r and r—i 2D colourmap with Pogson magnitude. The 
regression lines to radial colour {SBg — SBr) and (SBr — SBi), 



(SBg 

(SBr 



-SBr) 

' SB,) 



- + Clg- 

+ CIr-, 



■ log(a/a6o) 
log(a/a6o), 



(16) 



are calculated by linear least-squares fitting, and we define CIg-r 
and CIr-i as the radial colour gradient. The a is the equivalent dis- 
tance from galaxy centre derived by the axis ratio and position an- 
gle. The Petrosian 60% flux radius corresponds to an effective ra- 
dius with a pure de Vaucouleurs' profile. Generally the effective ra- 
dius is used for the normalization of the radial colour gradient, but 
profile fitting is needed to calculate the effective radius. However, 
some E+A morphologies are somewhat disturbed and it is difficult 
to obtain a stable profile fit. Thus, we use the Petrosian 60% flux 
semi-major axis for stability of calculation. When this linear least- 
squares fitting is applied, the a/aao of data points are restricted 
to 0.35 < a/aao < 1.0 to avoid a cusp in the galaxy centre and 
lower S/N regions in the outskirts. The minimum size of aeo in our 
samples applied to the colour gradient analysis is 2.83"(7.14 pix- 
els), and a/aao — 0.35 corresponds to the diameter of 2.0", since 
the typical F WHM of the SDSS point-spread function is ~ 1.5" 
iShimasaktTet aH2 00D. We confirmed by our tests using a control 
sample that the effects of seeing on the derived colour gradients 
are not significant. The (FWHM/2)/a6o v.s. |C/g_,.| shows nei- 
ther correlation (only 0.17 of correlation coefficient) nor systematic 
effect of seeing. 



4 RESULTS 

4.1 E-i-A Morphologies 

Since E-l-A galaxies have experienced starburst truncation fairly 
recently (< IGyr), E-l-A galaxies might still hold some traces in 
their morphology (e.g., dynamically disturbed signs). Therefore we 
might obtain some hint of the origins of E-l-A galaxies by examin- 
ing their morphology. 

First, we can confirm that E-l-A galaxies are centrally concen- 
trated using concentration index Ce - In Figure|2| we show the dis- 
tribution of Ce in 22 E-l-A galaxies. Galaxies that follow de Vau- 
couleurs' law give Ce = 0.29 and those with the exponential pro- 
file give 0.44, so Ce of early-type and late-type galaxies are clus- 
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Figure 3. g-, r- and i-band images, g—r and r—i 2D colourmaps (eacit top) and the rest colour profiles (eadi bottom) of our all E+A galaxies. The eiTor bars 
represent standard deviation of g—r or r—i values on sampling points. The details of explanation of image panels are inlaid in Figure^ The panels are placed 
in ascending order with respect to D4ooo- We numbered our E+As in this order, so this Figure includes No. l(top left), No.2(top right), No.3(foof/om left) and 
T^oA(hottom right) E+A galaxies. E+As with smaller D4000 tend to show positive slopes of radial colour gradients (bluer gradients toward the centre) and 
blue irregular patterns in 2D colourmaps. Meanwhile, those with larger Z)4ooo show negative or flat slope and moderate property in 2D colourmap. 
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Figure 3 - continued No.5, No.6, No.7 and No.8 E+A galaxies 
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Figure 3 - continued No. 21 and No.22 E+A galaxies 
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Figure 2. The distributions of concentration index Ce of 22 E+A galax- 
ies. Galaxies that follow de Vaucouleurs' law give Ce = 0.29 and those 
with the exponential profile give 0.44. All our E+A galaxies except one are 
centrally concentrated galaxies. 



tered on each point, respectively (See Yamauchi et al. 2005). Ce 
of 21 E+A galaxies are distributed in the same way as early-type 
galaxies, a nd the median Ce = 0.33 is a typical valu e of an E or 
SO galaxy iShimasaku et alJ200ll : IStrateva et al.l200lh . Our visual 
inspection^ of the SDSS images in TableQalso shows that our E+A 
galaxies are predominantly early-type galaxies. We can interpret 
these results as a suggestion that all except one of our E+A galax- 
ies are bulge-dominated systems. 

We showed one of the most splendid SDSS E+A galaxies with 
a tidal feature in Figure □ The SDSS imaging does not have high 
resolution (typical seeing size is ~1.5") compared with the HST or 

^ One of us (C.Y.) has performed a visual classification of the E+A 
galaxies in reference to the Third Reference Catalogue of Bright Galaxies 



| ;alaxies in reference to the 1 hird 1 
de Vaucouleurs et alll99il RC3). 



8m-class telescopes, but our 22 nearby or large E+A galaxies make 
a show of dramatic tidal tails or conspicuous disturbed morpholo- 
gies. We present all r-band negative images of our E+A galax- 
ies in Figure |3| on a logarithmic flux scale with 0% to 50%. The 
images in Figure |3| are placed in ascending order with respect to 
4000A break(i340oo), and we numbered our E+A galaxies accord- 
ing to this order. The physical size of aeo is displayed on the i- 
band image (right of the negative r-band image). To examine E+A 
morphologies, we present an image of a normal elliptical galaxy in 
Figure|4] Comparison between Figures|4]and|3|helps our examina- 
tion of E+A morphologies. These 22 E+A galaxies exhibit a variety 
of disturbed features, ranging from what could visually be classi- 
fied as a normal elliptical galaxy without disturbance (E+A No. 7) 
to one with impressive tidal features (E+A No.l). We notice that 
at least half of our E+A galaxies have traces of merger/interaction, 
and exhibit tidal features or disturbed morphologies. If we observe 
our E+A galaxies with a higher resolution using an 8m-class tele- 
scope, we should be able to examine the details of these features 
and the fraction of E+A galaxies which leave traces of dynamically 
disturbed signs might increase. 



4.2 2D Colour Properties: Radial Colour Gradients and 
Colour Morphologies 

An example of g—r and r—i 2D colourmaps and their radial colour 
profiles of normal early-type galaxy are presented in Figure|4| The 
3 positive images in the top left panel show g-, r- and i-band images 
by counts with seeing correction and resampling, and the 4 images 
on the right are g—r and r~i 2D colourmaps in asinh magnitude. 
The two rightmost panels of 2D colourmaps are pixel-to-pixel K- 
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Figure 4. Example of g-, r- and i-band image, g—r and r—i 2D 
colourmaps (top) and the rest colour profiles (bottom) of a normal early-type 
galaxy. 2D colourmaps are based on asinh magnitude, but Pogson magni- 
tude is used for colour profiles. The ii'-corrected 2D colourmaps are on the 
right, another has no ii'-correction applied. The error bars represent stan- 
dard deviation of g—r or r—i values on sampling points within an annulus. 
The details of explanation are inlaid in the top panel. The 2D colourmaps of 
ElHptical galaxy are featureless, and many of them have moderate negative 
slope (redder gradients toward the centre) at g—r and quite flat r—i colour 
gradient. 




Figure 5. The histograms of g—r radial colour gradient, showing E-l-A (top) 
and early-type (bottom) galaxies. 
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Figure 6. The histograms of r—i radial colour gradient, showing E+A (top) 
and early-type (bottom) galaxies. 



corrected, while the others are not. The seeing sizes of all images 
are indicated in the g-band image. The bottom panel shows g—r 
and r—i radial colour profiles using Pogson magnitude, and seeing 
size(FWHM/2) and physical scale of ago are displayed on the up- 
per abscissa. In this figure, we inlay the explanation of images and 
2D colourmap, and the same form is used in Figure |3| where we 
present all of our E-l-A galaxies. 

The 2D colourmaps of elliptical galaxies are basically fea- 
tureless, and the major difference between if-corrected and un- 
corrected is small. The radial colour profiles derived from 2D rest 
colour information also show moderate properties. Many elliptical 
galaxies have quite flat r—i colour gradients and moderate nega- 
tive slopes at g—r. Although there are many kinds of g—r radial 
gradients in early-type galaxies, the positive slope (bluer gradients 
toward the centre) at g—r profile is seldom seen. 

Meanwhile, many of our E-l-A galaxies show mysterious 2D 
colourmaps and radial profiles. We present 2D colourmaps and ra- 
dial profiles of all our E-l-A galaxies in Figure |3| The images and 
colour profiles in Figure |3| are placed in ascending order with re- 



spect to the 4000A break(D4ooo)> and we have numbered our E-l-A 
galaxies according to this order. The slopes of g—r radial profile 
of a number of E+A galaxies are positive (bluer gradients toward 
the centre) which is seldom seen in normal early-type galaxies. We 
list our 22 E-l-A galaxies in TableQwith CIg-r and CIr-i results. It 
turns out that ^--^63% of our E-l-A galaxies have a positive slope of 
radial g—r colour gradient. In the case of the r—i colour gradient, 
~77% of E-l-As have a positive slope. In addition, some E-l-A galax- 
ies show an irregular pattern in their g—r and r—i 2D colourmap 
without i^-correction — 'Colour Morphology' — asymmetrical 
and clumpy bluer region. The No. 1, No. 2 and No.3 E-l-A galax- 
ies are especially conspicuous; the bluer regions are shown near 
the galaxy centre, but some of the bluest regions are placed slightly 
away from nucleus. However, the A'-corrected 2D colourmaps tend 
to weaken this asymmetrical feature. We expect that such a region 
has a particular SED, and this may be a limit of isT-correction by 
fitting template SEDs. These 'irregularities' on 2D colourmaps are 
observed in E-l-A galaxies with 'positive slopes'. But not all of the 
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E+A galaxies with positive slopes necessarily liave conspicuous ir- 
regularities on the 2D colourmap. That is, the E+A galaxies with 
'positive slopes' of radial colour gradient tend to show irregular 
patterns on the 2D colourmap. 

To compare the radial colour gradients of the E+A galaxies 
with those of normal early-type galaxies, the distributions of radial 
colour gradients CIg-r and CIr-i are displayed in Figures|5|and|6| 
It is clear that normal early-type galaxies have negative g—r colour 
gradients, and almost all are placed within —0.2 < CIg-r < 0, and 
r~i colour gradients are concentrated at CIr-i ~ 0. Each median is 
—0.10 and 0.005, respectively. On the other hand, CIg-r and CIr-i 
of E+A galaxies are dispersed, and galaxies with large positive 
slopes are obviously increased. Each median is also shifted to posi- 
tive, 0.14 and 0.03, respectively. We then applied a Kolomogorov- 
Smirnov (K-S) two-sample test to find the probability that the two 
samples are drawn from the same parent distribution. The results 
show that g—r and r—i colour gradient distributions of E+A are 
different from those of early-type galaxies with a more than 99.99% 
significance level. 



4.3 Comparison between Radial Colour Gradient and Other 
Photometric/Spectroscopic Properties 

We showed that E+A galaxies have interesting radial colour gra- 
dients. Comparison between radial colour gradient and other pho- 
tometric/spectroscopic properties may provide a clue to the E+A 
evolution scenario. 

Figures Q and |8| show the relation between radial colour gra- 
dients CI and rest g—r or r—i colours. Observed colours in this 
paper are shifted to the restframe using the /f-correction software 
(v3_2) by Blantonetal. (2003). Our 22 E+A and 49 early-type 
galaxies are plotted using asterisks and open lozenges, respectively. 
The error bars on the top right indicate a typical error. Coefficient 
P on the bottom left is the Spearman linear correlation coefficient, 
and the solid line represents a linear least-squares fit using E+A 
data points. A data point with an open circle is No.l E+A galaxy 
rejected for P and linear least-squares fit, since the errors of CIg-r 
and CIr-i are too large (0.67 and 0.20, respectively). We expect 
that E+A galaxy No. 1 has too heavy a disturbance feature to permit 
calculation of the appropriate Petrosian radius. We can find obvious 
correlations of P ~ —0.7 level between radial colour gradients and 
colours. But some of the derived regressions appear to rely heavily 
on the one point, g—r — 0.48 E+A galaxy. So we have exam- 
ined the Spearman linear correlation coefficient P with 0.5 < g—r 
points only in Figure Q and gotten P = —0.62 and —0.59, re- 
spectively. On the other hand, the early-type galaxies show only 
P ~ —0.32 at the maximum absolute value in these diagrams. 
Therefore, derived P of E+As which indicate correlations between 
colour gradients and colours have some dependence upon the one 
point; however, it is not a prime factor. The results also support the 
idea that the photometric properties of E+A are significantly differ- 
ent from those of early-type galaxies. 

Next, we compare the radial colour gradients with spectro- 
scopic properties. In Figure |9| we plot radial colour gradients 
against the 4000A break (-D4000) which is sensitive to old stellar 
populations. Symbols and rejections are the same as in previous 
figures. This diagram evidently indicates the correlation between 
radial colour gradients and the 4000A break. Early-type galaxies 
show less than \P\ ~ 0.2 level correlations like those of colour gra- 
dients and colours. But normal early-type galaxies are not placed 
under the regression line, different from plots of colour gradients 



E+A * 




g-i- 

Figure 7. The g—r radial colour gradient plotted against g—r rest colour 
(top) and r—i radial colour gradient plotted against g—r rest colour (hot- 
torn), showing 22 E+A and 49 early-type galaxies. E+A and early-type 
galaxies are plotted using asterisks and open lozenges, respectively. The co- 
efficient P on the bottom left is the Spearman lineal' congelation coefficient, 
and the solid line represents a linear least-squares fit using E+A data points. 
The circled asterisks are No.l E+A galaxy rejected for P and linear least- 
squai'es fit. Error bars at the top right are typical errors of the observational 
data. 



against colours. -D4000 is shown to separate the two populations 
more clearly than the broad-band colours. 

To examine the relation between radial colour gradients and 
the amount of young A-type stars, we plot radial colour gradients 
against H5 EW in Figure [Tol Symbols are the same as in previous 
figures. Compared with the case of D4000. we cannot find a tight 
correlation. However, these panels show a trend for E+A galaxies 
with large positive slope of radial colour gradient to show strong 
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r—i 

Figure 8. The g—r radial colour gradient plotted against r—i rest colour 
(top) and r—i radial colour gradient plotted against r—i rest colour (bot- 
tom), showing 22 E+A and 49 early-type galaxies. E+A and early-type 
galaxies are plotted using asterisks and open lozenges, respectively. The co- 
efficient P at the bottom left is the Spearman linear coirelation coefficient, 
and the solid line represents a linear least-squares fit using E-l-A data points. 
The circled asterisks are No. 1 E-l-A galaxy rejected for P and linear least- 
squares fit. EiTor bars at the top right are typical errors of the observational 
data. 

H(5 EWs. However, this result might be affected by the cut-off of 
a/aeo < 0.35 for measuring radial colour gradient. 

4.4 Evolution Scenario for E-i-A Radial Colour Gradients 

In this subsection, we compare SED models and observational 
quantities to find out the evolution sce nario for E-l-A radial colour 
gradients. Using the GISSEL model bv lBruzual & Chariot 1,2003), 
we simulated th ree representative star formation histories follow- 
ing [Oot^ |20o3): (i) The Burst model, which has an instantaneous 




Figure 9. The g—r radial colour gradient plotted against D4000 (top) and 
r—i radial colour gradient plotted against D4000 (bottom), showing 22 
E-hA and 49 early-type galaxies. E-l-A and early-type galaxies are plotted 
using asterisks and open lozenges, respectively. Coefficient P on the bot- 
toms left is the Spearman linear correlation coefficient, and the solid line 
represents a linear least-squares fit using E-hA data points. The circled as- 
terisks are No. 1 E-l-A galaxy rejected for P and linear least-squares fit. Error 
bars at the top right are typical errors of the observational data. 



starburst (with a duration of 1 Gyr) at the beginning and no star for- 
mation thereafter, (ii) Constant star formation, (iii) Exponentially 
decaying star formation (with r = 1 G yr). In all three models, we 
use the Salpeter initial mass function i lSalDetej|ll9.5.5lb . Figure fTTI 
plots H(5 EWs against time (or galaxy age) for the three models. 
The dot-dashed, solid and dotted lines show the models with instan- 
taneous burst, constant star formation and exponentially decaying 
star formation rate. The burst model has a strong HJ EW right af- 
ter its truncation at 1 Gyr. However, its H(5 EW declines rapidly, 
and becomes less than 3A at 1 Gyr after the truncation. The expo- 
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Figure 10. The g—r radial colour gradient plotted against H(5 EWs {top) 
and r—i radial colour gradient plotted against H(5 EWs {bottom), showing 
22 E+A galaxies. The circled asterisks are No. 1 E+A galaxy. Error bars at 
the bottom right are typical eiTors of the observational data. 



nentially decaying model maintains a strong H(5 for a longer time. 
Its H(5 EW becomes 3A in 5 Gyrs. The constantly star-forming 
model maintains a large H(5 EW (>6A) beyond 13 Gyr. Although 
shown here for a comparison purpose, we omit the constantly star- 
forming model from our main discussion, since our E+A galaxies 
are selected to have no on-going star formation activity. 

We compare the models and the data on the H(5 EW v.s. D4000 
plane in Figure fT^ A caveat, however, is that models become less 
accurate on the plane since both H(5 EW and D4000 are more dif- 
ficult quantities to reproduce than broad band colours. For models, 
we use H(5 EWs given in the GISSEL model, which were mea- 
sured using the flux between 4083.50 and 4122.25A. This IS es- 
sentially the same window as used to measure H5 EWs from the 
observational data between 4082 to 4122A. For D4000, the SED 
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Figure 11. H<5 EWs are plotted against time (age) for three star formation 
histories with the GISSEL model. The dot-dashed, solid and dotted lines 
show the models with instantaneous burst, constant star formation and ex- 
ponentially decaying star formation rate. The models in this figure assume 
Salpeter IMF and solar metallicity. The inlaid panel is an enlarged plot of 
the burst model with 3 < H(5 EW < 9. 
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Figure 12. H<5 EWs are plotted against D4000 for the models with Z = 
0.0001(0.5% solar) and Z = 0.02(solar). Star formation histories are the 
burst and exponentially decreasing, shown by the dot-dashed and dotted, 
respectively. Observational data are plotted using asterisks. The circled as- 
terisk is No. 1 E+A galaxy. Error bars at the top right are typical errors of 
the observ ational data. T he contours show the distribution of all DRl 94770 
galaxies in lGotdEooS) . 
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Figure 13. The u—g colour is plotted against H5 EWs for the models 
with Z = 0.0001(0.5% solar) and Z = 0.02(solar). The dot-dashed and 
dotted lines are for the models with instantaneous burst and exponentially 
decaying star formation rate. Observational data are plotted using asterisks. 
The circled asterisk is No. 1 E+A galaxy. Error bars at the top right are 
typical errors of the observat ional data. T he contours show the distribution 
of all DRl 94770 galaxies in lGotd I2OO3I) . 



other hand, the u—g colour is preferable for examination of on- 
going star formation since it shows better agreement between the 
SED models and the observation data. In Figure fTsl we plot u—g 
against H(5 EW for the two models and the observational data. 
Symbols are the same as previous figures. The model u—g colour 
depends on metallicity to some extent, showing bluer u—g colour 
with decreasing metallicity. Our E-l-A galaxies are distributed be- 
tween the instantaneous burst model with Z = 0.001 and that with 
solar metallicity. Neither of the exponentially decaying models can 
explain the observed data. This indicates that E-l-A galaxies are in a 
post-starburst phase, and cannot be explained by more normal star 
formation histories. 

If the burst model is assumed, our young E-l-A galaxies with 
H(5 EW ~ SA evolve into those with RS EW ~ SA within 
~300 Myr by Figure fTTI Our radial colour gradient correlates with 
D4000 in Figure|9| and E-l-A galaxies with large positive slopes of 
radial colour gradient tend to show strong US EWs (Figure fTol . 
Our E-l-A galaxies are placed under a single track on the D4000 - 
H(5 EW plane; therefore, we can interpret these results as some 
E-l-A metamorphosis of large positive slopes into flat or negative 
slopes of radial colour gradients within the time scale. Taken to- 
gether, our observational results suggest that E-l-A galaxies have 
positive colour gradients, large US EW and small D4000 at the be- 
ginning, then evolve into negative colour gradients, small US EW 
and large D4000. It is revealing that the morphological quantity, 
colour gradients, correlates well with spectral properties such as 
RS EWand D4000. 



models uses the flux ratio of the 3750-3950A window to the 4050- 
5250A window (Bruzual 1983). Observationally, D4000 is mea- 
sured using the ratio of the flux in the 375 1 -3951 A window to that 
in the 405 1-4251 A jStoughton et alj|2002l) . We regard these two 
£^4000 measurements as essentially the same. 

In Figure fT2l the dot-dashed and dotted lines are for the mod- 
els with the instantaneous burst and the exponentially decaying 
star formation rate. We subtracted lA from the model H(5 EWs 
to compensate for possible stellar absorption. Different line widths 
are f or two different metallicities {Z = 0.001 and 0.02). lOotd 
{2003 ) showed that the models might have a slight shift toward 
larger D4000 and H(5 directions, compared with the distribution of 
all the observed galaxies. In the models, emission fillings are not 
included, and might be a factor in gaps between observational data 
and the model. However, the behavior of the models on this plane 
well reproduce expected behavior of galaxies, in the sense that 
star-forming galaxies evolve into large 1)4000 and small US EWs. 
Therefore, we regard the qualitative interpretation based on the 
models as valid. Observed data of E-l-A galaxies are plotted using 
asterisks, and we find for the most part that our E-l-A galaxies are 
located under a single evolution track. But it is uncertain whether 
these E-l-A galaxies run along the track of instantaneous burst or 
not, since models might have a slight shift reported by Goto (2003). 

iGotd i2003h also reported that it is more difficult to repro- 
duce emission lines with the current version of the model. The 
model gives the amount of Lya photons, so we could use it to es- 
timate the amount of Balmer or [Oil] line emissions. However, we 
found that the uncertainty in this process was too large to produce 
a realistic value comparable with the observational values (We do 
not know the escaping factor of Lya photons, etc.). We found a 
qualitative agreement in the long scale, but in the liner scale, the 
model [Oil] EW (or Ha EW) was far from agreement. On the 



5 DISCUSSION 

Based on SDSS broad-band imaging, we investigated spatial prop- 
erties such as morphologies, 2D colour properties and radial colour 
gradients of E-l-A galaxies. In this section, we discuss physical im- 
plications and the relevance of our results to other E-l-A-related 
studies. 

The concentration index, Ce, of our E-l-A galaxies was dis- 
tributed like early-type galaxies, and our visual inspection of the 
SDSS images also showed that our E-l-A galaxies are predom- 
inantly early-type galaxies. We can interpret these results as a 
suggestion that all of our E-l-A galaxies are bulge-dominated sys- 
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ous samples were contaminated by Ha-emitting galaxies. Perhaps 
since we selected 22 E-l-A galaxies with strict spectroscopic cri- 
teria, 'without Ha nor [Oil] emission', our 22 E-l-A galaxies are 
all bulge-dominated. Our results are also consistent with those by 
Ouintero (2004) and Blake et al. (2004). 

These findings have important implications for the physical 
origin of E-l-A galaxies when compared with theoretical simula- 
tions. It has long been known that an elliptical galaxy can be a final 
product of major merging (e.g., Barnes & Hemquist 1992), quite 
possibly accompanied by a tidal tai l feature in its early stage. In a 
more detailed major merger model. iBekkil il998t) showed that an- 
cestral equal-mass galaxies are completely destroyed so as to form 
a spheroid that looks more like an elliptical galaxy, and merging of 
two unequal-mass spirals is finally transformed into one SO galaxy 
that has a flattened oblate spheroid. We noticed that at least half of 
our E-l-A galaxies show tidal features or disturbed morphologies. 
Thus not only disturbed morphologies but also E/SO morphologies 
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Table 1. The list of E+A galaxies. The numbers are labeled in ascending order with respect to D4ooo- The digits in Name represent R.A. and Dec, for 
example, R.A. and Dec. of Name SDSSJ010617.76+140354.00(No.l6) are 01:06:17.76 and +14:03:54.00, respectively. The unit of aeo and line EWs are 
arcsec and A, respectively. The morphological type T = 0, 1, 2, 3, 4, 5, 6 and —1 represent E,SO,Sa,Sb,Sc,Sdm,Im and unclassified, respectively. 



No. 



Name 



2 



Mr u-g g-r r-i aeo Dtooo H(5 EW Ha EW [Oil] EW Ce T 



CIr 



16 SDSSJ010617. 

7 SDSSJ035652. 

18 SDSSJ083415. 

3 SDSSJ084918. 

19 SDSSJ093842. 

12 SDSSJ100743. 
6 SDSSJl 11050. 

4 SDSSJl 11 108. 
11 SDSSJl 15837. 

2 SDSSJ120418. 

13 SDSSJ120523. 
22 SDSSJ122256. 

20 SDSSJ132315. 

17 SDSSJ13335a 
9 SDSSJ140801. 

5 SDSSJ141003. 

21 SDSSJ141419. 
1 SDSSJ16133a 

15 SDSSJ162702. 
10 SDSSJ164608. 

8 SDSSJ170636. 

14 SDSSJ211348. 



76+140354.00 
44-061031.22 
60+375157.24 
96+462252.68 
96+000148.68 
68+554934.68 
88+005530.98 
16+004048.66 
68-021710.97 
96-001855.84 
28+643029.52 
16-011248.06 
60+630726.40 
64-001617.71 
68+514225.92 
60+603229.76 
20-031111.51 
24+510335.64 
64+432833.96 
16+363705.16 
48+334720.40 
24-063059.65 



0.038 
0.037 
0.168 
0.041 
0.091 
0.045 
0.152 
0.184 
0.088 
0.094 
0.082 
0.146 
0.175 
0.176 
0.160 
0.171 
0.047 
0.034 
0.046 
0.137 
0.124 
0.159 



■19.76 
■19.85 
■23.10 
■20.88 
■21.84 
■20.88 
■22.22 
■22.32 
■22.11 
■22.74 
■22.78 
■21.74 
■21.94 
■22.35 
■22.34 
■22.76 
■21.37 
■20.06 
■21.75 
■22.45 
■22.06 
■22.11 



1.33 0.67 

1.03 0.48 

1.23 0.89 

1.23 0.62 

1.40 0.76 

1.27 0.64 

1.21 0.81 
0.89 0.68 
1.27 0.69 
1.14 0.62 
1.23 0.66 
1.09 0.72 
1.35 0.97 
1.18 0.88 
1.17 0.83 
1.07 0.78 
1.60 0.79 
0.96 0.44 
1.54 0.77 
1.09 0.70 

1.22 0.74 
1.20 0.84 



0.34 2.96 
0.28 3.43 
0.41 4.93 
0.33 3.00 
0.38 4.07 
0.34 3.44 
0.39 3.34 
0.36 2.95 
0.36 3.35 
0.35 4.10 
0.36 4.64 
0.37 4.69 
0.43 3.00 
0.41 2.94 
0.39 2.88 
0.39 2.83 
0.38 5.53 
0.25 9.98 
0.37 5.38 
0.36 4.15 
0.37 3.08 
0.40 2.90 



1.51 5.53±0.32 
1.35 7.34±0.41 
1.51 5.58±0.23 

1.29 8.34±0.27 

1.51 5.57±0.47 

1.44 5.94±0.32 
1.33 6.34±0.32 

1.30 7.50±0.47 
1.43 6.10±0.28 
1.27 7.98±0.27 

1.45 6.79±0.29 
1.67 5.88±1.42 

1.52 6.90±0.93 
1.51 5.61±0.49 

1.38 6.81±0.37 
1.32 5.64±0.25 

1.53 5.81±0.34 
1.09 7.57±0.29 
1.50 5.69±0.22 

1.39 6.49±0.46 
1.37 5.99±0.39 
1.47 5.64±0.51 



2.10±0.08 
2.39±0.10 
0.74±0.02 
1.91 ±0.05 
1.1 5 ±0.04 
1.57±0.05 
1.28±0.05 
1.60±0.08 
1.75±0.07 
2. 11 ±0.06 
1.78±0.05 
1.19±0.05 
1.11 ±0.07 
2.21±0.14 
■0.12±-0.00 
0.67±0.05 
0.33±0.01 
0.54±0.02 
0.67±0.01 
1.12±0.03 
0.34±0.02 
0.96±0.03 



■0.13±-0. 
■0.46±-0. 
■0.22±-0. 
■1.17±-0. 
■1.62±-0. 
■2.06±-0. 
■0.94±-0. 
■0.57±-0. 

0.69±0. 
■1.18±-0. 

0.9 1±0. 

0.61±0. 

1.56±0. 

0.7 1±0. 

0.51±0. 
■1.08±-0. 
■l.ll±-0. 

1.10±0. 
■1.41±-0. 
■0.35±-0. 
■0.84±-0. 
■2.32±-0. 



02 0.28 
05 0.32 
,03 0.32 
,08 0.31 
21 0.36 
,20 0.28 
13 0.33 
08 0.36 
,09 0.34 
07 0.31 
,09 0.31 
,40 0.45 
,36 0.35 
15 0.31 
,07 0.33 
07 0.31 
18 0.34 
,21 0.36 
12 0.32 
04 0.35 
10 0.32 
56 0.33 



0.35±0.09 
0.64±0.21 
■0.14±0.34 
O.35±0.11 
■0.16±0.19 
0.35±0.12 
■0.00±0.18 
0.31±0.53 
0.16±0.09 
0.35±0.14 
0.30±0.16 
■0.32±0.45 
■O.30±0.33 
0.13±0.27 
0.34±0.23 
O.33±0.21 
■0.24±0.15 
■0.17±0.67 
0.08±0.12 
0.12±0.24 
0.20±0.23 
■0.06±0.18 



0.12±0.05 
0.27±0.09 
0.02±0.05 
0.12±0.03 
0.02±0.05 
O.08±O.O3 
0.03±0.04 
0.10±0.06 
0.03±0.01 
0.08±0.03 
0.06±0.03 
-0.12±0.09 
-0.04±0.05 
0.02±0.04 
0.09±0.03 
0.07±0.03 
-0.07±0.05 
-0.00±0.20 
0.02±0.03 
0.02±0.04 
0.05±0.04 
-0.03±0.03 



of our unalloyed E+A sample indicate possibilities of a merger ori- 
gin of E+A galaxies, and we can expedite the following discussion 
smoothly. 

We found irregular 'Colour Morphologies' - asymmetrical 
and clumpy patterns - at the centre of the g—r and r— i 2D 
colourmaps of E+A galaxies without if-correction. The HST ob- 
servations of_the five bluest E+A galaxies with z ~ 0.1 by 
lYaneetalJiiooi ) show the presence of asymmetric components in 
the residual images obtained by subtraction of the smooth and sym- 
metric model images from the data. The physical scale of the asym- 
metric components of the HST observation is comparable to that of 
our asymmetrical/clumpy pattern in the 2D colourmap. We expect 
that these irregularities of residual images and our 2D colourmap 
basically have the same origin, and one possible interpretation 
is that they are fluctuations created during the merger/interaction 
event, which also caused the central post-starburst phase. 

Our results showed that more than half of our E+A galax- 
ies exhibit the positive radial g—r and r—i colour gradients. The 
results of g—r radial colour gradients are in good agreement 
with a long-slit spe ctroscopic observation of 21 E+A galaxies by 
iNorton et alj i200lh who found that young stellar populations of 
E+A galaxies are more centrally concentrated than the older popu- 
lations. Although Bartholomew. Rose & Gaba (2001) reported that 
E+A galaxies on average tend to have slightly bluer radial gradi- 
ents toward the centre than the normal early-type galaxies, these 
E+A galaxies are detected in rich clusters and their spectroscopic 
criteria (3 < {US'yP EW> and -5 < [Oil] EW) are less strict 
than ours. It is perhaps difficult to compare our E+A galaxies 
with those selected with different criteria in different environments. 
However, we can find the same tendency for a substantial number 
of E+A galaxies to have a positive slope of radial colour gradient 
and for our strict spectroscopic criteria to produce clearer results 



of a K-S test (with a more than 99.99% significance level). The- 
oretically, models of interacting gas-rich spirals suggest that the 
most common result of such tidal perturbations is indeed to drive 
a large fraction of interstellar material close to the centre of each 
galaxy; the shocked gas rob bed of i ts angu lar momentum flows 
inward iSarnes & HernauisJ Il992h . iMihos . Richstone & Bothui] 
('1992^ showed not only that the vast majority of the star forma- 
tion arises in central regions of the merging galaxies, but also that 
flybys which involv e prograde disks cau s e sign ificant centralized 
starburst. Especiallv.l Mihos & Heanguistl i 19941) numerically sim- 
ulated galaxy mergers and predicted colour gradients in such sys- 
tems. In their Fig.7, they show that the colour gradient of merger 
remnants changes the slope of colour gradient from positive to neg- 
ative during 1-5 Gyr after the burst in a way qualitatively similar to 
our Fig|3| The agreement adds more credibility to the notion of the 
merger/interaction origin of E+A galaxies, and further justifies our 
interpretation of change of the slope as an evolutionary sequence. 
In summary, our findings of positive g—r radial colour gradients of 
E+A galaxies are consistent with the merger/interaction origin of 
E+A galaxies. 

We found a tight correlation between radial colour gradients 
and 4000A break. We also found that E+A galaxies follow a sin- 
gle evolution track on the -D4000 - H5 EW plane (in Fieure [T2l . 
Taken all together, these results can be interpreted in an evolution- 
ary sequence where E+A galaxies change their slope from posi- 
tive to negative, and they also change irregular into moderate of 
2D colourmap during the time scale of ~300 Myr. It is revealing 
that E+A galaxies show morphological metamorphosis (change in 
colour gradients and 2D colourmap) synchronously as their stellar 
population ages as indicated by 1)4000 and HJ EW. This may be 
the general picture of how E+A galaxies evolve after the trunca- 
tion of starburst possibly caused by the merger/interaction. It is a 
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concern that our radial colour gradients do not show a tight corre- 
lation with H5 EW, but the tendency of E+A galaxies with large 
positive slopes of radial colour gradient to show strong H5 EWs. 
However, this may be affected by our cut-off of the inner region 
(a/aeo < 0.35) to avoid seeing effects when measuring radial 
colour gradients. Observation with high-resolution imaging might 
remove these effects, and produce a more accurate picture. 



6 CONCLUSIONS 

We have investigated morphologies, 2D colour property and radial 
colour gradients of 22 true E+A galaxies with 5.5A < H5 EW < 
8.5A using the SDSS DR2 imaging data, in order to reveal evo- 
lution of E-l-A galaxies in terms of internal galaxy structures. We 
present this summary of our results: 

• Concentration index, Ce, and our visual inspection of the 
SDSS images suggested that our E-l-A galaxies are predominantly 
bulge-dominated systems, and at least half of our E-l-A galaxies ex- 
hibit tidal features or disturbed morphologies. 

• We found irregular 'Colour Morphologies' - asymmetrical 
and clumpy patterns - at the centre of g~~r and r—i 2D colourmaps 
of the E-hA galaxies. 

• We found that a substantial number of E-l-A galaxies have pos- 
itive slopes (bluer towards the centre) of radial g—r and r—i colour 
gradients. 

• Kolomogorov-Smimov two-sample tests showed that g—r 
and r—i colour gradient distributions of E-l-A galaxies are different 
from those of early-type galaxies with a more than 99.99% signifi- 
cance level. 

• Our E-l-A sample showed tight correlation between radial 
colour gradients and colours, and between radial colour gradients 
and 4000A break. We also found the tendency for E-l-A galaxies 
with large positive slopes of radial colour gradient to show strong 

EWs. 

• The compar i son between the GISSEL model 
iBruzual & Charlol l2003h and E-l-A's observational quanti- 
ties, HJ EW, -D4000 and u—g colour, indicated that almost all our 
E-l-A galaxies are located along a single evolution track. Therefore, 
these results are interpreted as indicating E-l-A radial colour 
gradients evolve from positive into negative slopes, coherently 
with RS EW (8-5A) within ~300 Myr. 

Our findings support the notion that E-l-A galaxies are post- 
starburst galaxies caused by merger/interaction, having undergone 
a centralized starburst. In addition, we presented a constraint on 
E-l-A's spatial evolution properties with a physical time scale. De- 
tailed comparison with future theoretical simulations would ad- 
vance our knowledge of the origin of E-l-A galaxies further. 
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